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Abstract 

The particle number projected BCS (PBCS) approximation is tested against the exact solution 
of the SO (5) Richardson-Gaudin model for isovector pairing in a system of non-degenerate single 
particle orbits. Two isovector PBCS wave functions are considered. One is constructed as a 
single proton-neutron pair condensate, while the other corresponds to a product of a neutron pair 
condensate and a proton pair condensate. The PBCS equations are solved using a recurrence 
method and the analysis is performed for systems with an equal number of neutrons and protons 
distributed in a sequence of equally spaced 4- fold (spin-isospin) degenerate levels. The results 
show that although PBCS improves significantly over BCS, the agreement of PBCS with the exact 
solution is less satisfactory than in the case of the SU(2) Richardson model for pairing between 
like particles. 
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I. INTRODUCTION 

Neutron-proton (np) pairing is a longstanding issue in nuclear structure Despite 
many efforts, the specific fingerprints of these correlations in existing nuclear data are not 
yet clear, nor the appropriate theoretical tools for their correct treatment. For many years 
the theoretical framework commonly used to describe the np pairing correlations was the 
generalized HFB approach [2]. In this approach the np pairing, both isovector and isoscalar, 
is treated simultaneously with neutron-neutron (nn) and proton-proton (pp) pairing. How- 
ever, although the generalized BCS approach treats on equal footing all type of pairing 
correlations, most of BCS calculations show that they rarely mix Thus, in general, there 
are three BCS solutions which seem to exclude each other: one with nn and pp pairs; the 
second, degenerate to the first in even-even N = Z nuclei, with isovector np pairs; and the 
third with isoscalar np pairs. 

Various studies have shown that the restoration of particle and isospin symmetries and the 
inclusion of higher order correlations improve significantly the predictions of BCS approach 

nnnpfi 

for systems with np pairing |5|, Ig, |7|, |8|. To restore exactly these symmetries, projection 
operators or projected generator coordinate methods are commonly employed 9|. Less 
discussed in the literature is an alternative method based on the recurrence relations satisfied 
by the isovector pairing Hamiltonian averaged on projected BCS (PBCS) wave functions. 
In this paper we will implement this method to analyze the dependence of isovector pairing 
correlations on particle number conservation. As trial wave functions we will use two PBCS 
condensates, one formed by isovector np pairs and another by nn and pp pairs. Contrary 
to the BCS approximation for a system of an even number of pairs, the PBCS solutions 
corresponding to these two pair condensates are not degenerate. To analyze how much 
these PBCS solutions could improve over the generalized BCS approach will shall use the 
exactly solvable SO (5) Richardson-Gaudin pairing model |10| . Several previous studies have 
been carried out in the one-level degenerate S0(5) model These studies clarified the 
limitations of the BCS approximation, and the corresponding extensions taking into account 
pair fluctuations in the RPA formalism or using boson expansion theories [s], Q, [l2|. Studies 
on number and isospin projection on the isovector pairing Hamiltonian with non-degenerate 
single-particle levels have been reported in ^]. However, these studies were tested against 
a solution proposed by Richardson [isl and later on shown to be incorrect for systems 
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with more than two pairs [l^. The exact solution of the non-degenerate isovector pairing 
Hamiltonian has been given by Links et al. [15] and afterwards generahzed to seniority 
non-zero states, arbitrary degeneracies, and symmetry breaking Hamiltonians in lOj]. This 
solution will be used here as a benchmark to test the accuracy of PBCS approximations for 
describing the isovector pairing correlations. 



II. FORMALISM 



We will consider an isovector (T = 1) pairing Hamiltonian with a constant pairing 
strength 

^ = E - 9 E V(J. + 1/2)(J.' + 1/2)P,+ P,,., (1) 

irriT i,i',T 

where = ':^['^X^'j'^ol is the isovector pair creation operator. The first column in the 
couplings refers to total angular momentum and the second column to total isospin. 

The Hamiltonian (1) is a particular example of the exactly solvable S0(5) Richardson- 
Gaudin integrable models. It is is exactly solvable for arbitrary single particle energies ej. 
and pair degeneracies ji + 1/2. The exact solution of these class of Hamiltonians has been 
given in Ref. [l^. Here we will treat a simplified version for a system of L equidistant single- 
particle levels of pair degeneracy 1, that is jj = 1/2. The exact solution for this system will 
be used as a test for the PBCS approximation with isovector pairing. For comparison we 
shall also show the results of the proton-neutron BCS approximation. The generalized BCS 
model used in this paper is described in Ref. |l^. As in the case of a single degenerate 
level Q], within the BCS approximation the Hamiltonian (1) has two solutions: (A) a BCS 
solution with a non-zero proton-neutron gap, A„p ^ 0, and zero gaps for neutron-neutron 
and proton-proton pairs, i.e., A„ = Ap = 0; (5 ) a BCS solution with A„ = Ap 7^ and 
Ap„ = 0. The two solutions (^4 ) and {B ) exclude each other and are degenerate in energy 
for a system with an even number of pairs. In the next section we will present the PBCS 
equations corresponding to these two BCS solutions. The PBCS formalism will be given 
in the form of recurrence relations, and it can be applied to general (density-independent) 
isovector pairing interactions, irrespectively of whether they are integrable or not. 
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A. PBCS approximation with isovector proton-neutron pairs 



We shall first consider a PBCS wave function corresponding to the solution (A), i.e., 
formed by N isovector neutron-proton pairs. It has the following form 

\^>-W.^^of\0>, (2) 
where Fq is the collective neutron-proton pair operator 

L 

^+ = J2^^P^l (3) 

1=1 

This wave function is not normalized and the factor in front is chosen to simplify the form 
of PBCS equations. The mixing amplitudes Xi are determined by minimizing the energy 
functional 

^ <N\H\N> 

^ ' <N\N> ^ ' 

The norm and the expectation vahic of the Hamiltonian are calculated by using recurrence 

relations. Thus, it can be shown that the norm of the wave function (2) satisfies the equation 
<7V|7V>=l^a;^<Ar-l|Ar-l>--L^^3<7V-l|i^+|7V-2> (5) 

i 



where 



< A^|i^^ |iV - 1 >= a;^ < AT - l|iV - 1 > -L^ < TV - 1\P,^\N - 2 > . (6) 

To get the norm corresponding to the system with N proton-neutron pairs the equations 

above should be iterated starting with < 1|1 >= < l|-Pjo|0 

The expectation values of the particle number operators Ni, which give the occupation 
probabilities of the single-particle levels, can be calculated from the equation 

< N\Ni\N >= 2xi < N\P+\N - 1 > (7) 

where the matrix elements in the r.h.s. are given by Eq.6. 

Finally, the matrix elements of the pairing force are given by the equations 

< N\P^^Pjo\N > = ^xy^ <N- 2\P;,P,o\N - 2 > 

+Xj < N\P+\N - 1 > -^x^jXi <N- 1\P+q\N - 2 > 

+Sij^\< N -2\N -2> -l < N - 2\Ni\N - 2 >] 
■^4 2 



< N\P+P,, + Pl,P,-i\N > = ^ < iV - 2|i^|P,i + ,P._i|iV - 2 > 

x4 1 
+5ij^[< A^-2||A^-2 > — <N -2\Ni\N -2 >] 

These equations above are iterated starting from < l\P^PjQ\l >— XiXj and < l|Pj^Pji + 
P,t,Pj-i\l>^ 0. 

B. PBCS approximation with proton-proton and neutron-neutron pairs 

We will now consider a PBCS wave function corresponding to the BCS solution (B), i.e., 
given by a product of two condensates formed by nn and pp pairs. This trial wave function 
has the form 

\MM >^ \M > ®|M >= ^(r+r+)^|o > (8) 

where M denotes the number of nn and pp pairs, M — N/2, while r+ and r+ are the 
collective pair operators for neutrons and protons (see Eq.(9) below). As defined here, the 
wave function (8) is well suited for cvcn-cven nuclei. For odd-odd nuclei the corresponding 
wave function is formed by M — {N — l)/2 neutron-neutron and proton-proton pairs plus 
two unpaired nucleons that block the corresponding levels affecting the pairing correlations. 

Since the Hamiltonian (1) is symmetric in isospin, ior N — Z systems the collective 
proton and neutron pair operators should have the same mixing amplitudes, i.e., 

L L 

r+ = 5]y,/^|, v;^Y.y^pU (9) 

Due to the same reason, the norms for the neutron and proton wave functions and the 
matrix elements for the neutron-neutron and proton-proton interaction should satisfy similar 
recurrence relations. Therefore below we shall give only the recurrence relations for one kind 

of particles, i.e., neutrons. Thus, the norm of the neutron state \M > and the average of 
neutron number are given by 

<M|M>=^5;]l/,^<M-l|M-l>-^ J]yf <M-l|i^+|M-2> (10) 

i i 

< M\Ni\M >^ 2yi < M\P+\M - I > (11) 

where 

< M\P+\M -l>^yi<M -l\M -l> -y1 <M - l\P^\M - 2 > . (12) 
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The matrix elements of the neutron-neutron pairing interaction are given by the equations 

< M\PlP,^\M > = yly] <M- 2|P+P,i|M - 2 > 

+yj < M\P+\M - 1 > -y^Vi <M- l|P+|Af - 2 > 
+6ijyt[< M - 2\M - 2 > - < M - 2\Ni\M - 2 >] (13) 

The iterations are started with the matrix elements < l|Pj^Pji|l >= yiyj. Eqs. (10-13) are 
very simple and can be used as an alternative to the projecting operator method commonly 



applied for systems with like-particle pairing 



17|. 



The matrix elements of the T = 1 proton-neutron interaction involve the total wave 
function \MM >. They are given by the recurrence relation 

< MM\P^^Pjo\MM >= yfy] < M - IM - l\P+Pjo\M - IM - 1 > + 
6ijxt <M - l\M -l>[< M -l\M -l> - < M - l\Ni\M - 1 >] 

The starting matrix elements are < ll|Pj^Pjo|ll >= 5ijX^. As can be seen from the equa- 
tions above, the recurrence relations for the PBCS wave functions (2) and (8) are very similar 
and easy to implement in numerical calculations. 



III. RESULTS AND DISCUSSIONS 

The results presented in this section correspond to a sequence of L equally spaced 4- 
fold degenerate levels (total angular momentum j = 1/2) with single particle energies Ei = 
[i — l)/2, z = 1,2, ...L and filled with = L/2 proton-neutron pairs (quarter filling). We 
have considered systems with A^ = 2 to A^ = 12 pairs, which correspond to typical sizes 
of open shell N = Z nuclei. The strength of the pairing interaction is varied to cover 
all regimes from weak to strong coupling. For these systems we will test the accuracy of 
the PBCS approximations comparing correlation energies, odd-even mass differences and 
occupation probabilities against the exact solution. We will start this comparison focusing 
on correlation energies. They are defined as 

Ecorrig) = Enorig) " E{g) (14) 

where Enor and E{g) are the ground state energies of the system in the normal and in 
the correlated phase respectively. Some representative results are shown in Figs. 1-3. All 
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FIG. 1: Correlation energy for 4 pn pairs 

energies are given in units of the single particle level spacing. In these figures PBCO cor- 
responds to the variational wave function (2) of = np pairs, and PBCSl corresponds 
to the variational wave function (8) of nn (T^ = l)and pp (Tj = —1) pairs. The two BCS 
solutions corresponding to these two types of pairs are called BCSO and BCSl. In even-even 
systems these two BCS solutions are degenerate and are called simply BCS. Particle number 
projection breaks this degeneracy. 

As can be see in Figs. 1 and 2, both PBCS solutions perform better than BCS for even 
systems, with PBCSl capturing more correlations and lowering the ground state energy. On 
the other hand, as shown in Fig. 3, for a system with an odd number of pairs the lowest 
energy solution is PBCSO. It can be also seen that due to the blocking, in the systems 
with odd number of pairs the solution PBCSl becomes higher in energy even than the BCS 
solution. 

The errors relative to the exact results are shown in Fig. 4. It can be seen that although 
PBCS gives better results than BCS, the errors remains significant. The reason is that the 
PBCS functions (2) and (8) do not take into account properly the pairing interaction among 
the pairs with a different from what is considered in the trial wave function. For example, 
let's consider the systems with 8 and 7 pairs and the interaction strength g = 0.4. In the 
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FIG. 2: Correlation energy for 8 pn pairs 

system with 8 pairs the wave function PBCSl gives an energy of -25.71 for the = ±1 
part of the hamiltonian compared to -0.84 for the = part. The situation is opposite for 
the system with 7 pairs: in this case the wave function PBCSO gives an interaction energy 
of -18.04 for the = component compared to about -1.42 for = ±1. 

Another quantity we have analyzed is the odd-even mass difference along the N = Z line 
defined as 

A^^\M) = ^[2E{M + 1) - E{M) - E{M + 2)]. (15) 

Fig. 5 shows the odd-even mass difference for a system with M = 8 pn pairs as a function of 
interaction strength. It can be seen that the PBCS results start to deviate significantly from 
the exact values when the interaction becomes stronger. How the odd-even mass difference 
depends on the number of pairs is depicted in Fig. 6. As expected, the BCS results do not 
show the staggering exhibited by the exact solution. This is because in BCS the solutions 
(A) and (B) are degenerate in energy. On the other hand the staggering is present in the 
PBCS calculations. This is due to the fact that going from the even-even to odd-odd systems 
the ground state is changing from PBCSO to PBCl, which are not degenerate. As seen in 
Fig. 6, the shift between the two solutions overestimates the oscillations present in the exact 
solution. The reason is that the errors in odd systems are larger than in even systems (see 
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FIG. 3: Correlation energy for 7 pn pairs 

Fig. 4). 

Next we shall discuss shortly the occupation probabilities corresponding to BCS and 
PBCS calculations. Figs. 7-8 show the quantity = vf{l — vf), where vf is the occupation 
probability of the orbit i. In BCS is the pairing tensor and determines the pair transfer 
form factor. From Figs. 7-8 we can see that PBCS gives results close to the exact solution 
for both values of the coupling strength. BCS overestimates the value of at the weak 
coupling (g=0.25) in the region around the Fermin energy, where the pairing correlations are 
stronger. Conversely, the states further than an energy interval of the order of the pairing 
gap are underestimated. These results are similar to the ones obtained in Ref. [iSj] for 
like-particle pairing. For stronger interactions (g=0.4) BCS gives results closer to the exact 
solution. 

Up to now we have considered two distinct PBCS wave functions. The question is if one 
could get extra binding by mixing together the wave functions PBCSO and PBCSl. This 
is indeed what happens for a system formed by two pn pairs. The results are shown in 
the table below. As can be seen, by mixing the two PBCS states one gets practically the 
exact result for the correlation energy. As expected, we get the extra binding when the 
mixed state has the total isospin equal to zero. For systems with more than two pairs a trial 
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FIG. 4: Errors of the correlation energies in systems with 7 and 8 pn pairs 

wave function with zero isospin cannot be constructed by mixing only the states PBCSO 
and PBCSl. Consequently for such systems we do not get extra binding by mixing the two 
PBCS wave functions. 

TABLE I: Correlation energies for a system composed of two isovector pn pairs distributed in four 
levels with the energies ei = {i — l)/2, i = 1,2, 3, 4. The binding energies and the pairing strength 
g are in units equal to the distance between two consecutive single-particle levels. 
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Exact 


PbcsO+Pbcsl 


Pbcsl 


PbcsO 


0.1 


0.05587 


0.0557 


0.0376 


0.0189 


0.2 


0.22006 


0.2192 


0.1517 


0.0779 


0.4 


0.81330 


0.8114 


0.5924 


0.3233 


0.6 


1.64761 


1.6461 


1.2551 


0.7364 


0.8 


2.61989 


2.6190 


2.0601 


1.2972 


1.0 


3.66946 


3.6689 


2.9487 


1.9683 
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FIG. 5: Odd-even mass difference for a system with 8 pn pairs 
IV. SUMMARY AND CONCLUSIONS 

We have analyzed the accuracy of PBCS approximation for describing isovector pairing 
correlations in N = Z systems. The study was done for an exactly solvable hamiltonian with 
S0(5) symmetry. In the PBCS calculations we considered two kind of trial wave functions: 
(1) a condensate of isovector neutron-proton pairs; (2) a product of two condensates formed 
by neutron-neutron and proton-proton pairs. The solution (1) gives the lowest ground state 
energy for odd-odd N = Z systems while the solution (2) provides the lowest energy for 
even-even systems. The PBCS approximation gives much better correlation energies than 
BCS, and it is able to describe the staggering of odd-even mass difference calculated along 
the N = Z line. However, compared to the pairing between like particles, for which the 



PBCS approximation give results very close to the exact solution of the SU(2) model jl8| . 
the accuracy of PBCS approximation for isovector pairing is less satisfactory. The reason is 
that the PBCS is not able to treat correctly that part of the isovector force which describes 
the interaction among the pairs which are not included in the PBCS condensate. Going 
beyond PBCS would imply including the isospin projection and/or taking into account 
quartet correlations. We are currently working along the later direction. 
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FIG. 6: Odd-even mass difference calculated along N=Z line 
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FIG. 7: = vf{l — vf ) as a function of the single-particle energies £i for a system with 8 pn pairs 
and pairing strength g=0.25 
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FIG. 8: The same as in Fig. 8 for the pairing strength g=0.4 
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